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A Cascadic Multigrid Method for GPE Problem* 


Xiaole Han^ and Hehu Xie^ 


Abstract 

A cascadic multigrid method is proposed for the GPE problem based on the 
multilevel correction scheme. With this new scheme, the ground state eigen¬ 
value problem on the finest space can be solved by smoothing steps on a series 
of multilevel finite element spaces and some nonlinear eigenvalue problem solv¬ 
ing on a very low-dimensional space. Choosing the appropriate sequence of 
finite element spaces and the number of smoothing steps, the optimal con¬ 
vergence rate with the optimal computational work can be arrived. Some 
numerical experiments are presented to validate our theoretical analysis. 

Keywords. Bose-Einstein condensation; Gross-Pitaevskii equation; mul¬ 
tilevel correction; cascadic multigrid; nonlinear eigenvalue problem; finite el¬ 
ement method. 
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1 Introduction 

The aim of this paper is to design a cascadic type multigrid finite element method for 
solving Gross-Pitaevskii equation (GPE) which is a time independent Schrodinger 
equation. The hnite element method for GPE problem and general semilinear eigen¬ 
value problem has been investigated by [3, 15]. The corresponding error estimates 
are also given. 

Recently, a type of multilevel correction method is proposed to solve eigenvalue 
problems in [6, 12, 13]. In this multilevel correction scheme, the solution of eigen¬ 
value problem on the hnal level mesh can be reduced to a series of solutions of 
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boundary value problems on the multilevel meshes and a series of solutions of the 
eigenvalue problem in a very low-dimensional space. Therefore, the cost of com¬ 
putation work can be reduced largely. A multigrid method for the GPE has been 
proposed in [14] where a superapproximate property is founded. Therefore, the aim 
of this paper is to construct a cascadic multigrid method to solve the ground state 
solution of Bose-Eienstein condensates (BEG). The cascadic multigrid method for 
second order elliptic eigenvalue problem is given in [5]. This method transforms the 
eigenvalue problem solving to a series of smoothing iteration steps on the sequence 
of meshes and eigenvalue problem solving on the coarsest mesh by the multilevel 
correction method. Similarly to the cascadic multigrid for the boundary value prob¬ 
lem [2, 9], we only do the smoothing steps for the involved boundary value problems 
by using the previous eigenpair approximation as the start value and the numbers 
of smoothing iteration steps need to be increased in the coarse levels. The order of 
the algebraic error for the final eigenpair approximation can arrive the same as the 
discretization error of the finite element method by organizing suitable numbers of 
smoothing iteration steps in different levels. The nonlinear eigenvalue problems on 
a very low-dimensional space are solved by self-consistent iteration or Newton type 
iteration which reduces the nonlinear eigenvalue problem to a series of linear ones. 

The rest of this paper is organized as follows. In the next section, we introduce 
the hnite element method for the ground state solution of BEG. A cascadic multigrid 
method for solving the non-dimensionalized GPE is presented and analyzed in Sec¬ 
tion 3. In Section 4, some numerical tests are presented to validate our theoretical 
analysis. Some concluding remarks are given in the last section. 


2 Finite element method for GPE problem 


This section is devoted to introducing some notation and the finite element method 
for the GPE problem. In this paper, we shall use the standard notation for Sobolev 
spaces and their associated norms and semi-norms (cf. [1]). For p = 2, we 

denote and HliVt) = {n G H^iVL) : v\gQ = 0 }, where = 0 is 

in the sense of trace, || • = || • ||s, 2 ,r 2 - The letter C (with or without subscripts) 

denotes a generic positive constant which may be different at its different occurrences 
through the paper. 

For simplicity, we consider the following non-dimensionalized GPE problem: Find 
(A, u) such that 


-Am -|- Wu -|- CkPw 


u 


In 


Am, in G, 

0, on 5G, 

1 , 


( 2 . 1 ) 


where G C TZ'^ {d = 2, 3) is a bounded domain with Lipschitz boundary dfl, ( is 
some positive constant and W(x) = 'jixl • -f > 0 with 71 , • • • , 7 d > 0 . 
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In order to use the finite element method to solve the eigenvalue problem (2.1), 
we need to define the corresponding variational form as follows: Find (A, m) eTZxV 
such that b{u, u) = 1 and 


a{u,v) = Xb{u,v), Vn G V, 


( 2 . 2 ) 


where V := Hq{Q) and 

a{u,v) := / [VuVv+ Wuv + (\u\‘^uv)d^l, b{u,v) := / uvdVt. (2.3) 

The existence, uniqueness and simplicity of the smallest eigenpair of eigenvalue 
problem (2.2) have been given in [3]. 

To simplify the notation, we also define H^iVL) inner-product a(-, •) as 

a{w,v) := / 'Vw'Vvd^l, \/w,v ^V. (2.4) 

Jn 

Now, let us define the finite element approximations of the problem (2.2). First we 
generate a shape-regular decomposition of the computing domain C 7^'^ (d = 2, 3) 
into triangles or rectangles for d = 2 (tetrahedrons or hexahedrons for d = 3). The 
diameter of a cell K E Th is denoted by hx and the mesh size h describes the 
maximum diameter of all cells K G Th- Based on the mesh Th, we can construct a 
finite element space denoted by 14 C V. For simplicity, we set 14 as the linear finite 
element space which is defined as follows 

Vh = {vh G 4(fi) I v^Ik G Vu WK g %}, (2.5) 

where Vi denotes the linear function space. 

The standard finite element scheme for eigenvalue problem (2.2) is: Find {Xh, Uh) G 
TZ X Vh such that b{uh, Uh) = 1 and 


a{uh,Vh) = Xhb{uh,Vh), WvhEVh. ( 2 . 6 ) 

Then we define 

4(m) := inf \\u-Vh\\i. (2.7) 

Vh&Vh 

Lemma 2.1. ([3, Theorem 1]) There exists h^ > 0, such that for all 0 < h < ho, 
the smallest eigenpair approximation {Xh,Uh) of (2.6) having the following error 
estimates 


u-Uh\\i < C6h{u), (2.8) 

U-Uh\\o < Cr]a{Vh)\\u-Uh\\i <CT]aiVh)5hiu), (2.9) 

|A-Aft| < C{\\u-Uh\\l + \\u-Uh\\o) <CT]aiVh)6h{u), (2.10) 
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where rjaiVh) is defined as follows 


Va{Vh) = sup inf ||r/-u,,||i (2.11) 

with the operator T being defined as follows: Find T f & such that 

a{Tf,v) + 2{C\u\^{Tf),v) - (A(T/),u) = {f,v), Wv e 

where = {v & Hlifil) : \ J^uvdQ = 0}. Here we use the fact that 5h{u) < 
CVaiVh). 


3 Cascadic multigrid method for GPE 

Recently, a multilevel correction scheme is introduced in [6, 12, 13] for solving 
Laplace eigenvalue problems. Based on their involved idea, we propose a type of 
cascadic multigrid method for GPE problem (2.2) in this paper. The main idea 
in this method is to approximate the underlying boundary value problems on each 
level by some simple smoothing iteration steps. In order to describe the cascadic 
multigrid method, we hrst introduce the sequence of hnite element spaces and the 
properties of the concerned smoothers. 

In order to design multigrid scheme, we hrst generate a coarse mesh Th with the 
mesh size H and the coarse linear hnite element space Vh is dehned on it. Then we 
dehne a sequence of triangulations Th^, of C determined as follows. Suppose 
Thi (produced from Th by regular rehnements) is given and let Th^. be obtained from 
Thf._^ via one regular rehnement step (produce subelements) such that 

hk ~ (3.1) 

where the positive number fi denotes the rehnement index and larger than 1 (always 
equals 2). Based on this sequence of meshes, we construct the corresponding nested 
linear hnite element spaces such that 

VffCK, CR,, (3.2) 

The sequence of hnite element spaces I4j C Vh 2 C • • • C Vh„ and the hnite element 
space Vh have the following relations of approximation accuracy 

ha(Rf/) > 4 i(m), dhfiu) fc = 2,---,n. (3.3) 

In fact, since the ground eigenvalue A of (2.2) is simple (see [3]) and the computing 
domain is convex, we have the following estimates 

ha(Rff) ~ H, rjaiVhfi) ^ hk and dhfiu) hk, k = ,n. (3.4) 
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Remark 3.1. The relation (3.3) is reasonable since we can choose Shf.{u) = hk {k = 
1, • • • ,n). Always the upper hound of the estimate Sh^{u) < hk holds. Recently, we 
also obtain the lower bound result Sh^{u) > hk (c.f. [7]). 

For generality, we introduce a smoothing operator Sh ■ Vh ^ 14 which satishes 
the following estimates 


r \\s(fwh\\i < ^j^\\wh\\o, 

{ ii^r^.iii < ikhiii, (3.5) 

I < \\ShWh\\i + wsjfvhWi, 

where 4 is a constant independent of h and a is some positive number depending 
on the choice of smoother. It is proved in [4, 8, 11] that the symmetric Gauss-Seidel, 
the SSOR, the damped Jacobi and the Richardson iteration are smoothers in the 
sense of (3.5) with parameter a = 1/2 and the conjugate-gradient iteration is the 
smoother with a = 1 (cf. [9, 10]). 

Then we dehne the following notation 


Wh = Smooth{Vh, /, fh, m, Sh) (3.6) 

as the smoothing process for the following boundary value problem 

a{uh,Vh) = b{f,Vh), yvh e 14, (3.7) 

where fh denote the initial value of the smoothing process, Sh denote the chosen 
smoothing operator, m the number of the iteration steps and Wh is the output of 
the smoothing process. 

Now, we come to introduce the cascadic multigrid method for the eigenvalue 
problem ( 2 . 2 ). Assume we have obtained an eigenpair approximations (A^'“,m^'') G 
n X 14^,. We design the following cascadic type one correction step to improve the 
accuracy of the current eigenpair approximation G 77 x Vh,. 

Algorithm 3.1. Cascadic type of One Correction Step 

1. Define the following auxiliary source problem: Find g such that 

+ 6 14.^.. (3.8) 

Perform the smoothing process (3.6) to obtain a new eigenfuction approxima¬ 
tion G 14^_^^ by 

= SmoothiVh,^,, - {W + C\u^'^\‘^)u’^Cu’^Cmk+i,Sh,^,). (3.9) 
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2. Define a new finite element space = Vh + span{M^''+i} and solve the 

following eigenvalue problem: Find (A^''+i, eIZx such that 

b^y^hk+i^ y^hk+i'^ = 1 and 

a{u^>‘+\v'^^+^) = e'l4'=+\ (3.10) 

Summarize the above two steps by defining 

(A^''+i,m^''+i) = SmoothCorrection{VH,Vh^^^, 

Based on the above algorithm, i.e., the cascadic type of one correction step, we 
can construct a cascadic multigrid method for GPE as follows: 

Algorithm 3.2. GPE Cascadic Multigrid Method 

1. Solve the following GPE problem in the initial finite element space Vh^: Find 
(A^bn^i) ^ TZ X 14j such that 

a{u’^\Vhfi) = \’^^h{u^\vhfi), Vufti e 14i. 

2. For k = 1, - ■ ■ ,77,-1, do the following iteration 

= SmoothCorrection{VH,Vhi^^-^, ,u^'",mk+i, Shk_^fi. 

End Do 

Finally, we obtain an eigenpair approximation (A^",m^") eTZx Vh„. 

In order to analyze the convergence of Algorithm 3.2, we introduce an auxiliary 
algorithm and then show its superapproximate property. Similarly, assume we have 
obtained an eigenpair approximations {Xh^,Uh^.) eTZx 14^.. The following auxiliary 
one correction step is dehned as follows. 

Algorithm 3.3. Auxiliary One Correction Step 

1. Solve the following auxiliary source problem: Find G such that 

= Xhf,b{uhf,-iVh fe + l ) 

-((fT + C|Mh,r)Mh„Vh,+J, e 14,^,. (3.11) 

2. Define a new finite element space VH,hk+i = '^h + span{M/i^_^4 + span{M^'=+i} 
and solve the following eigenvalue problem: Find (A/i^_^j ,E 71 x Vn^hk+i 
such that h{uh,,^^,Uhk+i) = 1 o,nd 

a{uh^,^^,VH,hk+f) =\t,+fi){uh^+i,VH,hk+i), e 147,/xfc+i- (3.12) 
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Summarize the above two steps by defining 

= AuxiliaryCorrection{VH, Xh^,Uhk,u’^’°^fi. 

Algorithm 3.4. GPE Auxiliary Multilevel Correction Method 

1. Solve the following GPE problem in the initial finite element space Vh-^: Eind 
{Xhi,Uh^) eTZx 14 j such that 


a{uh^,Vhfi) = Xhfi){uh^,Vhfi), Vuft, e 14i. 

2. Eor k = 1, - ■ ■ ,77,-1, do the following iteration 

= AuxiliaryCorrectioniVn, Xh^,Uh;^,u’^'‘+^). 


End Do 

Einally, we obtain an eigenpair approximation {Xh„,Uh„) eTZx Vh„. 

Before analyzing the convergence of Algorithm 3.2, we show a superapproximate 
property of Uh^ obtained by Algorithm 3.4. The similar result is also analyzed in 
[14]. 

Theorem 3.1. Assume Uh^ (k = 1, ■ ■ ■ ,n) are obtained by Algorithm S.f and Uhi^ 
(k = 1, ••• ,77^ the standard finite element solution in Vh^. If the seguence of fi¬ 
nite element spaces 14^, • • • ) ^hn the coarse finite element space Vh satisfy the 
following condition 

Cpa{VH)fi^ < 1, (3.13) 

the following estimate holds 

\\uh,,-Uhji < Cr]a{Vh^,)5hfiu), /c = l,---,77, (3.14) 

and 

\\uh,,-Uhjo < Cr]a{Vh^,)5hfiu), /c = l,---,77, (3.15) 

where C is a constant only depending on the eigenvalue X. The eigenvalue approxi¬ 
mations Xh^. and Xh^. have the following estimates 

\Xhu-'^hu\<Cria{Vhfi)5hfiu), fc = 1,2, • • • ,77. (3.16) 

Proof Dehne := \Xh^^ - Xhfi + \\uh^ - UhJlo, k = 1,2, - ■ ■ ,n. And it is obvious 
that Ehi = 0. From (2.6) and (3.11), we have 


O‘i'^hk+1 '^hk+i y '^hk+i ) 

= b{Xh,._^_iUhf._^_i — Xhf.Uhi^,Vhf,^fi 
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fc+1 ^hk'^hf, lloll'^^h fc + 1 111 

+ 11 ^/ 1 ^.^^ — 'W/ife||o(||'W/ife+i|lo,6,rj llo,6,f^) ll'*^^fe+i II l) 

< C{\\Xhk+iUhk+i - XhkUhkWo + \\XhkUhk - -^/ifeW/iJIo 

+ {\\Uh,+, -MhJ|o+ ll^h, -MhJ|o)(||Mh,+J|?+ ||M/.J|?))||Wh,+il|l 
^ C{jXhl,^l — Xh^. I + llw/l;.,,.! — Uhk ||o + ll'^/lfc + l II 1 - 

It leads to the following estimates 

||hhj._,_i — Uhf,j^^ II 1 — C'(|l/ifc+i ~ Xhf. I + llh/ij,^^ — Uhf, ||o + (3-17) 

Note that the eigenvalue problem (3.12) can be regarded as a finite dimensional 
subspace approximation of the eigenvalue problem (2.6). Similarly to Lemma 2.1 
(see [3]), from the second step in Algorithm 3.3, the following estimate holds 

llh/ij.,,.! —'U/ij,_,_i||i < inf — VH,hk+x\\i ^ (3.18) 

Then combining (3.17) and (3.18) leads to 

||'h/ii;+i ~ '^hk +1 II 1 — C'(l^/ife+i ~ Xhf, I + ll'h/ij.^j — Uhy. ||o + ^hi,) ■ (3.19) 

From the properties of 14^, C 14^.^^, Vn^hk C 14^., Lemma 2.1 and (3.3), we have 

\\uh^+^ -_UhJ\i < C6h^{u), \\uh^^^ - Uhjo < C'r/a(I4j||hft,+i - Uhji, 

\Xh^^-Xh^\ < Cr]a{Vh^)Sh^{u), \\uh^-UhJ\o<Cr]a{VH)\\uh^-UhJ\i, 
\Xhk-Xhk\ < Cr]a(yH)\\uh;^-UhJli- 

Substituting above inequalities into (3.19) leads to the following estimates 

||hh,+,< C{ria{Vh,)Sh,{u)+eh,) 

< C{r]aiVh^)6h^{u) + VaiVH)\\uh^ - UhJli). (3.20) 

When fc = 1, since Uh^ := Uh^ and Xh^ := Xh^, we have 

\\uh2 -Uh^Wi < Cr]a{Vh^)5hM)- ( 3 - 21 ) 

Based on (3.3), (3.20), (3.21) and recursive argument, we have the following esti¬ 
mates: 

k 

IIS*.-S/.JI 1 < (Cr,.)4,_,(«) 

i=2 

k 


< cf{Y, {Cva(v„)P^)'‘-’)v.(VkMu) 

i=2 

(JR‘2 

Therefore, the desired result (3.14) holds under the condition Cr]a{VH)l3‘^ < 1- Fur¬ 
thermore, (3.15) and (3.16) can be obtained directly from Lemma 2.1 and the prop¬ 
erty Vhm +1 c 14,^1. □ 


Note that C VH,hk, we can obtain the following estimates which play an 
important role in our analysis. 

Lemma 3.1. [3, Theorem 1] Let and Uh^, Vn^hk defined in Algorithms 

3.1 and 3.3. Then the following estimates hold: 

\\u^^-uhji < C\\uh,-u^^\\u (3.23) 

< C?7a(Vff)||M^'' - M/iJll, (3.24) 

< Cva{VH)\\u^^ -Uh,\\i. (3.25) 

Proof Since l4'= C Vhm^ according to (3.10) and (3.12), can be viewed as the 
spectral projection of Uhk- Then from Lemma 2.1 and the dehnitions of Vn^hk and 
1/^'“, we have 

\\uh, - < C inf 114/,,-n^'=||i<C' inf 

< C\\uh,-u^'^\\i, (3.26) 

which is the desired result (3.23). 

Similarly, we also have (3.24) by the following argument 

W^hk - < Cr|a{VH'^)\\Uh^, - u'^^Wi < C7]a{VH)\\Uhk “ 

Furthermore, (3.25) can be obtained directly from Lemma 2.1 and the proof is 
complete. □ 

Remark 3.2. Since Vh C and Vh C VH,hk, from Lemma 2.1, we have 

- MhJIi < 14^'“ - w||i + ||w - w/ij|i < C5h{u). (3.27) 


Now, we come to give error estimates for Algorithm 3.2. 

Theorem 3.2. Assume the eigenpair approximation is obtained by Algo¬ 

rithm 3.2, {XhnjUhfr is obtained by Algorithm 3.4 and the smoother selected in each 
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level Vhf. satisfy the smoothing property (3.5) for k = 1, - ■ ■ ,n. Under the conditions 
of Theorem 3.1, we have the following estimate 

(3.28) 

k=2 

and the corresponding eigenvalue error estimate 

<Crb(y„)^ (1 + (3 29) 

k=2 

Proof. Define eh^, '■= — Uh^ for fc = 1, • • • , n. Then it is easy to see that Ch^ =0. 

From Lemma 3.1, the following ineqnalities hold 

l|e/xfc+illi = - Mftfe+illi < C'l|Mh,+i - 

< + (3.30) 

For the hrst term in (3.30), together with (3.8), (3.11), Lemma 3.1 and (3.27), we 
have 


- Xh.UhJo 

< C{\X^>^ -Xh,\ + \\u^>^ -Uh,\\o) 

< CpaiyH)\\u^*‘ -Uhjl 

= C'ha(^r/)l|c/ifel|i- (3.31) 

For the second term in (3.30), due to (3.5) and (3.31), the following estimates hold 

< iisry‘(s'“« -s-jii. + iicyK. - “'“)iii 

< iis:;:r(S'**" -s*„.)iii+iisr“‘(“'..« - -“'“iii 

C 1 

< -l|w/xfc+i - W/^JIo + \\Uhu - 

^k+l ^fc+1 

< {l + Cpa{VH))\\ehu\\i+ -MfeJIo. (3.32) 

^k+i 

According to Lemma 2.1, (3.3), Theorem 3.1 and its proof, we have 


< Cpa{Vh^+f)dh^+Au)- (3.33) 

Combining (3.30), (3.31), (3.32), (3.33) and (3.4), we have 

\\eH+^\l<{^ + Cpa{VH))\\eh^\l + -^5h^^,{u), fc = l,---,n-l. (3.34) 


m 


fc+i 
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Based on (3.34), the fact e/j^ = 0 and the recursive argument, the following estimates 
hold 

lle^Jli < {l + C7]a{VH))\\eh^_A\i + 

< (1 + CT]aiVH)Y\\eh^_^\\i + (1 + Cr]a{VH))-^6h^_,{u) + -^ShM 

n ^ 

< c^(i + c,„(y„))"-‘—4.W. 

iff'].. 

k=2 ^ 


This is the desired result (3.28). The estimate (3.29) can be obtained from Lemma 
2.1 and (3.28). □ 


Corollary 3.1. Under the conditions of Theorem 3.2, we have the following esti¬ 
mates: 


n 

< C'(^a(14„)(5h„(M) + ^ 


{l + Cr^a{VH)) 


n—k 


k=2 

n 


mi 


ikjii) , (3.35) 




(1 + Ci).(Va)) 


n—k 


k=2 


mV 


iScJu) ■ (3.36) 


Now we come to estimate the computational work for Algorithm 3.2. Define the 
dimension of each linear finite element space as 


Nk := dim Vh^, k = !,■■■ ,n. 


Then we have 


Nk 




d{n—k) 


Nn, 


fc = 1, • • • ,n. 


(3.37) 


Different from the linear Laplace eigenvalue case, in the second step of Algorithm 
3.2, we have to solve a nonlinear eigenvalue problem on the newly constructed coarse 
space Vlf". Always, some type of nonlinear iteration method is used to solve this 
nonlinear eigenvalue problem. In each nonlinear iteration step, we need to assemble 
the stiff matrix on the finite element space {k = 2, • • • ,n), which needs the 
computational work 0{Nk). Fortunately, the matrix assembling can be carried out 
by the parallel way easily in the finite element space since it has no data transfer. 


From Theorem 3.2, in order to control the global error, it is required that the 
number of smoothing iterations in the coarser spaces should be larger than the fine 
spaces. To give a precise analysis for the final error and complexity estimates, we 
assume the following inequality holds for the number of smoothing iterations in each 
level mesh: 


hn / m" V / 


k = 2, - ■ ■ ,77,-1, 


(3.38) 
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where m = rrin, a > 1 and C > 1 are some appropriate constants. 

Now, we give the hnal error and the complexity estimates for Algorithm 3.2. 

Theorem 3.3. Under the conditions (3.3), (3.38) and + CH) < 1, for any 

given 7 G (0,1], the final error estimate 


\u 


'^hn 111 — 


holds if we take 


m > 



(3.39) 


(3.40) 


where Q = 1/(1 - ^^-^(l + CH)). 

Assume the GPE problem solved in the coarse spaces Vh and need work Mh 
and Mh^, respectively. We use P computing-nodes in Algorithm 3.2, and let w 
denote the nonlinear iteration times when we solve the nonlinear eigenvalue problem 
(3.10). IfC/<^ < d, the total computational work of Algorithm 3.2 can be bounded 

by o(^[l + j)Nn + M/ij + M/^log(A^„)j and furthermore 0{Nn) provided Mh 
N n, Mhi < Nn and ^ < C. While if C,/a = d, the total computational work 
can be bounded by 0{{1 + zu/p)Nnlog{Nn) + log(A'„)) and furthermore 

0{Nn log(A'„)) provided Mh N^, M^^ < and wjp < C. 


Proof. By Theorem 3.2, together with (3.1), (3.4), (3.28), (3.38) and (3^ ^(1+CiJ) < 
1 , we have the following estimates 




k=2 

n 




k=2 

n 


hr. 


< c ^(1 + 


k=2 


n—2 


< c 


cjfiEW-a+cH-)) 

k=0 

hn 1 


< 


1 - /3i-C(l + CH) 

CC,^ 


'-hr 




(3.41) 


Then it is obvious that we can obtain — Mft„||i < 'yhn when m satishes the 
condition (3.40). 
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Let W denote the whole computational work of Algorithm 3.2, Wk the work on 
the fc-th level for k = 1, ■ ■ ■ , n. Based on the dehnition of Algorithms 3.1 and 3.2, 
(3.1), (3.37) and (3.38), the following estimates hold 


W = E Wk E Mill A ^ ^ + 'y ^ loS/3(-^n 


k=l 


k=2 


k=2 


, , _ 1 / / 1 \ in-k){d-C/a) 

< Mhi + CMh log(A'„) + maA“A'„ ^ j 


k=2 


■W / 1 \ d{n-k) 


k=2 


/S 


< M,, + C'M^,log(iV„) + C'-iV„+maA“iV„5^M 

^ k=2 


1 \ in-k){d-(/a) 


Then we know that the computational work W can be bounded by 0(^Mhi + 
MH^og{Nn) + (l + 'aj/p)Nn) when d — (/a > 0 and by 0{Miii + MH^og{Nn) + 
(1 + tt7/p)A'n log(A^n)) when d — (/a = 0. It is also obvious they can be bounded by 
0{Nn) and C>(A'„ log(A'„)), respectively, if Mh A'n, Mh^ < Nn and wjp < C are 
provided. □ 

Remark 3.3. Since we have a good enough initial solution u^'=+i in the second step 
of Algorithm 3.2, solving the nonlinear eigenvalue problem (3.10) always dose not 
need many nonlinear iteration times (always zu < 3). 

Corollary 3.2. Under the same conditions of Theorem 3.3 and (d.fO), ifChn < 7, 
we have the following estimate 


- UhJ\i < ‘^ihn. (3.42) 

If we choose the conjugate gradient method as the smoothing operator, then 0 = 1 
and the computational work of Algorithm 3.2 can be bounded by 0{fl + w /p)Nn + 
Mhi + Alfj^og^Nn)) or 0{Nn) provided Mh Nn, M^i < Nn and wjp < C for 
both d = 2 and d = 3 when we choose 1 < C < d. 

When the symmetric Gauss-Seidel, the SSOR, the damped Jacobi or the Richard¬ 
son iteration acts as the smoothing operator, we know a = 1/2. Then the computa¬ 
tional work of Algorithm 3.2 can be bounded by 0{{l+vj /p)Nn+Miii+MH log(A'„)) 
{0{Nn) provided Mh <C Nn, Mf^ < Nn and wjp < C) only for d = 3 when we 
choose 1 < C < 3/2. In the case of a = 1/2 and d = 2, from Theorem 3.3 
and its proof, we can only choose (/ = 1 and then the hnal error has the esti¬ 
mate ||uA _ < Chn{l + computational work can only be 

bounded by O ((1 w/p)Nn log(A'„) Mhi + Mh log(A^„)) (C>(A^n log(A^n)) provided 
Mh < Nn, Mhi < Nn and wjp < C). 
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4 Numerical exapmle 


In this section, we give a numerical example to illustrate the efficiency of the cascadic 
multigrid scheme (Algorithm 3.2) proposed in this paper. Here, we choose the 
conjugate-gradient iteration as the smoothing operator (a = 1) and the number of 
iteration steps by 


ruk = for = 2, • • • , n 

with m = 2, a = 2, = 2, (C = 1.8 and [r] denoting the smallest integer which is 

not less than r 

Here we give the numerical results of the cascadic multigrid scheme for GPE 
problem on the two dimensional domain = (0,1) x (0,1) with W = x\ + x\ 
and C = 1- The sequence of finite element spaces are constructed by using linear 
element on the series of meshes which are produced by the regular refinement with 
/9 = 2 (connecting the midpoints of each edge). In this example, we use two meshes 
which are generated by Delaunay method as the initial mesh Th^ and set Th = Th-^ 
to investigate the convergence behaviors. Figure 1 shows the corresponding initial 
meshes: one is coarse and the other is fine. 

Algorithm 3.2 is applied to solve the GPE problem. For comparison, we also solve 
the GPE problem by the direct finite element method. From the error estimate result 
of GPEs by the finite element method, we have 

5h{u) K. h, r]a{Vh) ^ h. 

Then from Gorollary 3.2, the following estimates hold 

\\Uh„ - u’^'h < Chn, \\Uh^ - Wo < CHhn, < CHhn- 

We consider the Delaunay meshes (see Figure 1). 




Figure 1: The coarse and fine initial meshes for the unit square (left: H=l/6 and right: 
H=l/12) 
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Figure 2 gives the corresponding numerical results for the GPE problem on the 
initial mesh illustrated by the left mesh in Figure 1. The corresponding numerical 
results for the GPE problem on the initial mesh illustrated by the right mesh in 
Figure 1 are shown in Figure 3. 


Eigenvalue errors for H=1/6 Eigenfunction errors for H=1/6 



Figure 2: The errors of the cascadic multigrid algorithm for the GPE problem, where 
and denote the eigenfunction and eigenvalue approximation by Algorithm 3.2, and 
and denote the eigenfunction and eigenvalue approximation by direct eigenvalue 
solving (The left figure is the eigenvalue errors and the right figure is the eigenfunction 
errors which both correspond to the left mesh in Figure 1) 


Eigenvalue errors for H=1/12 Eigenfunction errors for H=1/12 




Figure 3: The errors of the cascadic multigrid algorithm for the GPE problem, where 
and A^ denote the eigenfunction and eigenvalue approximation by Algorithm 3.2, and 
and A^'"^ denote the eigenfunction and eigenvalue approximation by direct eigenvalue 
solving (The left figure is the eigenvalue errors and the right figure is the eigenfunction 
errors which both correspond to the right mesh in Figure 1) 

From Figures 2 and 3, we hud the cascadic multigrid scheme can obtain the same 
optimal error estimates as the direct eigenvalue solving method for the eigenfunction 
approximations in the hf^-norm. 

Remark 4.1. Note that by (3.36) and (3.42), we do not prove the optimal conver¬ 
gence rate for eigenvalue error (i.e. |A/i„ — A^"| < Ch(^). However, it is shown in 
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the left of Figures 2 and 3 that < Ch"^. 

5 Concluding remarks 

In this paper, we present a type of cascadic mnltigrid method for GPE problem based 
on the combination of the cascadic mnltigrid for boundary value problems and the 
multilevel correction scheme for eigenvalue problems. The optimality of the com¬ 
putational efficiency has been demonstrated by theoretical analysis and numerical 
examples. 
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